Surface anchoring on liquid crystalline polymer brushes 
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Abstract 

We present a Monte Carlo study of the surface anchoring of a nematic fluid on swollen layers of grafted liquid 
crystalline chain molecules. The liquid crystalline particles are modeled by soft repulsive ellipsoids, and the chains 
are made of the same particles. An appropriately modified version of the configurational bias Monte Carlo algorithm 
is introduced, which removes and redistributes chain bonds rather than whole monomers. With this algorithm, a 
wide range of grafting densities could be studied. 

The substrate is chosen such that it favors a planar orientation (parallel to the surface). Depending on the 
grafting density, we find three anchoring regimes: planar, tilted, and perpendicular alignment. At low grafting 
densities, the alignment is mainly driven by the substrate. At high grafting densities, the substrate gradually loses 
its influence and the alignment is determined by the structure of the interface between the brush and the pure 
solvent instead. 
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1. Introduction 

Liquid crystals have fascinated condensed mat- 
ter physicists, chemists and material scientists for 
many decades [1,2] - first, because they exhibit a 
large number of beautiful phases with intriguing 
symmetries and interesting material properties, 
and second, because of their technological use, 
particularly in the domain of display devices. One 
important aspect for technical applications is the 
so-called surface anchoring [3-5]: Surfaces orient 
nearby molecules, and these in turn orient the 
whole bulk of the liquid crystal. Alignment lay- 
ers are key components of liquid crystal display 



devices [5]. In particular, one is interested in de- 
signing surfaces which orient liquid crystals in a 
well-defined way in any desired direction. 

Possible candidates are liquid crystalline poly- 
mer brushes [6-10] . Halperin and Williams [6] have 
suggested to create a situation where a substrate 
which favors planar anchoring (parallel to the sur- 
face) competes with a stretched main-chain liquid- 
crystalline brush which favors homeotropic anchor- 
ing (perpendicular to the surface). They predicted 
that this conflict would result in tilted anchoring 
above a critical grafting density S c . In their sce- 
nario, the transition between planar and tilted an- 
choring is continuous. Above the transition, the 
grafting density can be used to tune the anchoring 
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angle. 

Experimentally, dense liquid crystalline polymer 
brushes have been prepared [9], and their align- 
ment properties have been studied [10]. So far, the 
conditions in the experiments were such that the 
substrate and the brush favor the same type of 
alignment. However, a competition as required by 
the scenario of Halperin and Williams can be in- 
troduced by a straightforward modification of the 
substrate [10]. 

Motivated by that work, we have investigated 
the consequences of such a competition by Monte 
Carlo simulations of an idealized liquid crystal 
model. We have studied brushes of liquid crys- 
talline polymers in a nematic solvent on a substrate 
which favors planar anchoring. The relaxation 
times in such systems are very high. Therefore, we 
have implemented a new version of the configura- 
tional bias Monte Carlo algorithm [11—14]: Chains 
are broken up such that their monomers turn into 
solvent particles, and solvent particles are linked 
together to form new chains. With this algorithm, 
a wide range of grafting densities could be con- 
sidered. We shall see that the resulting anchoring 
scenario turns out to be quite rich. 

We introduce the simulation model in the next 
section, and describe the algorithm in the third 
section. The results are presented and discussed in 
section 4. We summarize and conclude in the last 
section. 



2. The Simulation Model 

Our system consists of soft ellipsoidal particles 
with elongation k = cr cnd _ cni /<r sidc _ sidc = 3. Two 
particles i and j with orientations and Uj sep- 
arated by the center-center vector interact via 
the purely repulsive pair potential 
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with x = (« 2 — 1)/( k2 + !)• The function Uij ap- 
proximates the contact distance between the two 
ellipsoids in the direction = r^/r^ [15]. 

Solvent and chain particles have the same non- 
bonded interactions. In addition, chain monomers 
are connected by anharmonic springs of length b 
with equilibrium length bo and a logarithmic cut- 
off at \b — bo I = b s . At \b — b | < b s , the spring 
potential is given by 



T4(b) = -fc s /2b 2 ln[l-(b-b ) 2 /b 2 )]. 



(3) 



The direction of a bond is given by the unit vector 
b. A stiffness potential 

t4(u,bi,b 2 ) = -k a (|ubi| + |ub 2 | + 2bib 2 ) (4) 

is imposed, which penalizes non-zero angles be- 
tween the orientation u of a monomer and its two 
adjacent bonds t>i and b 2 . The system is confined 
by hard walls at two sides z = and z = L z : The 
distance between the walls and the centers of el- 
lipsoids, which have an angle 9 with respect to the 
surface normal, must exceed 



d z {6) = a /2 v/l + cos 2 ^) ( K 2 - 1). 



(5) 
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The chains are grafted to the walls at one end, the 
grafting points are arranged on a regular square 
lattice. 

We have studied systems with roughly N — 2000 
solvent particles (the numbers varied in the differ- 
ent runs) and up to 242 chains of four monomers. 
The simulation boxes were rectangular of size Lu x 
L|| x L z with periodic boundary conditions in the 
lateral directions and fixed boundary conditions 
in the z direction. The thickness L z was allowed 
to fluctuate, so that the simulations could be per- 
formed at constant pressure. The lateral size i|| = 
12<To was kept fixed in order to maintain a constant 
grafting density S. Trial moves included monomer 
displacements, monomer rotations, rescaling of L z , 
and the configurational bias moves described in 
the next section. The initial configuration was set 
up such that all particles pointed in the x direc- 
tion. The system was then equilibrated over at least 
1 million Monte Carlo steps, and data were col- 
lected over 5 or more million Monte Carlo steps. 
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One "Monte Carlo step" includes on average N at- 
tempts of monomer displacements, 2N attempts of 
monomer rotations, one attempt of L z rescaling, 
and one attempt of a configurational bias move. 

The model parameters chosen were k s — k a — 
10., 60 = 4. and b s = 0.8. Here and throughout, 
we use scaled units defined in terms of e , a , and 
the Boltzmann constant ks- Unless stated other- 
wise, the simulations were performed at constant 
temperature T = 0.5 and constant pressure P = 
3. The average number density in the bulk was 
(p) = 0.313. 



3. Configurational Bias Monte Carlo in a 
Solvent 

The most frequently used continuous space ver- 
sion of the configurational bias Monte Carlo al- 
gorithm [12,16] removes chains from the system 
and inserts new monomers into the system. How- 
ever, this method is not effective in our system, 
because the space is densely filled with monomer 
and solvent particles, and little free volume for 
new particles is available. Therefore, we have re- 
sorted to a different scheme. Instead of remov- 
ing the monomers of a chain, we only remove the 
connecting bonds and redistribute them such that 
new chains are formed from solvent particles. This 
method resembles closely the lattice version of the 
configurational bias Monte Carlo algorithm [12]. 
The only difference is that the lattice is now ran- 
dom, with "lattice sites" given by the centers of 
solvent or monomer particles. 

In practice, we proceed as follows: First, we re- 
move all the bonds from a randomly chosen chain. 
Then we find all m solvent particles, which are 
at a distance be [bo — b s , b a + b s ] from the graft- 
ing point, and choose one of them according to the 
Boltzmann probability 
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Here b\ is the length of the bond between the graft- 
ing point and the particle, bi its direction, and Ui 
the orientation of the particle. The normalization 
factor Z\ is the sum 
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over the Boltzmann factors for all n\ potential can- 
didates j . We continue in the same spirit and con- 
struct a new trial chain monomer by monomer: 
Given the (a — l)th monomer, we identify all n a 
solvent particles at the distance b G [bo — b s , bo + b s ] 
from the (a — l)th monomer, and pick one with 
the probability 

P a = ^exp{-^[V s (b a ) 

~ ^a(|u Q — ib Q I + 2b«_ib a + |u Q b Q |)]} (8) 
with 
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The total probability of constructing a given new 
chain with M monomers is 
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For the final Metropolis step, one also has to com- 
pute the probability P i d that the original, initially 
removed chain is constructed. The way to do this is 
exactly analogous, except that the monomers are 
not picked randomly. The new chain is accepted 
with the probability 

HW^ now = min (1, exp(-A£/T)), (11) 

-tnew 

where AE is the energy difference between the old 
and the new chain. Since no particles were moved, 
AE only contains contributions of the bond ener- 
gies. If the move is not accepted, the original chain 
is restored. 

A similar method has been described recently by 
Wijmans et al [14]. The implementation of this al- 
gorithm was essential for the success of our simula- 
tions. With the new algorithm, the configurations 
equilibrated within a million Monte Carlo steps 
and reliable data for a range of grafting densities 
could be collected. 
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Fig. 1. Bulk phase diagram in pressure-temperature space, 
calculated from simulations with three different system 
sizes as indicated. 
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Fig. 2. Solvent in contact with the bare substrate; config- 
uration snapshot (top) and number density profile p vs. z 
(bottom). 

4. Results 

Before investigating the full system, we have 
characterized the different constituents of our 
system separately. 

The bulk phase diagram of the solvent parti- 
cles was determined from simulations of cubic sys- 
tems with periodic boundary conditions in all di- 




Fig. 3. Configuration snapshot of grafted chains without 
solvent at the grafting density E = 0.34/o-y. 
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Fig. 4. Grafted chains without solvent: Average tilt an- 
gle (9) of the head-to-end vector of the chains as a function 
of the grafting density E. 

rections. The resulting phase diagram in P — T 
space is shown in Figure 1. The state point of our 
later simulations at T = 0.5, P = 3 is well in the 
nematic phase. The coexistence densities at P = 3 
are p No matic « 0.29, and /cWropic ~ 0.28. The sys- 
tems were too small to allow for an accurate de- 
termination of the coexistence gap in p — T space. 
In P — T space, however, the phase boundaries are 
not affected much by finite size effects (Figure 1). 

A bare substrate without grafted chains orients 
the particles parallel to itself. A snapshot is shown 
in Figure 2. The density profile reveals strong lay- 
ering effects in close vicinity to the wall, which de- 
cay rapidly in the bulk. 

Finally, we have studied the behavior of brushes 
without solvent. A snapshot for a system with the 
moderate grafting density £ = 0.34/cr§ is shown 
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in Figure 3. The chains are disordered, with chain 
orientations distributed around 9 — (6 being the 
angle with respect to the substrate normal). The 
average angle (8) of the head-to-end vectors de- 
creases rapidly with increasing grafting density £ 
(Figure 4). 

We turn to the discussion of the results for the 
full system. Figure 5 shows a snapshot of a brush 
with solvent at the grafting density £ = 0.34/ctq. 
Comparing this figure with Figure 4, we note that 
the presence of the solvent makes the brush much 
more compact, i.e., the solvent particles introduce 
an effective attractive interaction between the 
chains. The chains are tilted in a common direc- 
tion (Figure 5 (top)). The tilt propagates into the 
bulk of the system outside of the brush (Figure 
5 (bottom)). As anticipated in the introduction, 
we thus find that a brush can induce tilt in the 
orientation of a nematic solvent. 

The comparison of Figure 4 and Figure 6 (top) 
shows that the average tilt angle (6) of the head- 




Fig. 5. Configuration snapshot of grafted chains with sol- 
vent at grafting density E = 0.34/o-q. Top: Just chain 
monomers (solvent particles are transparent); Bottom: 
Chain monomers (light) and solvent particles (dark). 



to-end vector of chains increases in the presence of 
the solvent at all grafting densities. Most notably, 
the chains remained tilted up to the highest graft- 
ing densities that we considered. However, the tilt 
of the chains does not always produce tilted an- 
choring in the nematic bulk fluid. Figure 6 (bot- 
tom) demonstrates that the solvent particles have 
less tilt than the chain particles, and stand almost 
perpendicular to the surface at the highest graft- 
ing density. In fact, a closer inspection of our data 
shows that the remaining average angle (9) stems 
mainly from the contribution of the solvent parti- 
cles inside of the brush, and from orientational dis- 
order of the particles outside of the brush. A more 
detailed discussion will be presented elsewhere [17] . 

Figure 6 (bottom) reveals a second feature of the 
system: At low grafting densities, the wall aligns 
the solvent particles in a planar way, i. e., the 
alignment is basically determined by the anchoring 
properties of the bare substrate. At the grafting 
density £ « 0.12, the tilt angle of solvent particles 
jumps discontinuously to a lower value. The data 
thus indicate that there is a first order transition 
between planar and tilted anchoring. 



5. Discussion and Summary 

Our simulation data disclose the following sce- 
nario: We find three different anchoring regimes. 
At low grafting densities, the anchoring on the 
wall is planar. At high grafting densities, it is 
homeotropic, i.e., perpendicular to the wall. Tilted 
anchoring can be achieved in a range of interme- 
diate grafting densities. The transition between 
planar and tilted anchoring is presumably first 
order. The transition between tilted and perpen- 
dicular anchoring has not been studied in detail, 
because simulations at high grafting densities 
are very expensive. The data are consistent with 
theoretical mean field considerations [18], which 
suggest that it is continuous. According to theory, 
the anchoring behavior at high grafting densities 
is controlled by the structure of the interface be- 
tween the polymer layer and the nematic fluid. 
The competition between attractive chain interac- 
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tions, the translational entropy of the solvent, and 
its elasticity, drives a continuous transition from 
tilted to homeotropic alignment. 

From a practical point of view, low graft- 
ing densities are more interesting, because they 
can be realized more easily in experiments. The 
anchoring transition which we observe in that 
regime resembles that predicted by Halperin and 
Williams [6]. However, the underlying mechanism 
is probably quite different. The theory of Halperin 
and Williams relies on the existence of hairpins, 
abrupt reversals in the directions of the chains, 
whereas our "polymers" are much too short to 
support such defects. Nevertheless, the transition 
seems to be controlled by the organization of the 
chains in the polymer layer. We shall investigate 
this in more detail in a future publication [17]. 
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Fig. 6. Average tilt angle (9) as a function of the graft- 
ing density S. Top: Head-to-tail vector of grafted chains; 
Bottom: Single particles. 
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